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Abstract 

The fusion probability in "hot" fusion reactions leading to the synthesis of super-heavy nuclei is 
investigated systematically. The quasi-fission barrier influences the formation of the super-heavy 
nucleus around the "island of stability" in addition to the shell correction. Based on the quasi- 
fission barrier height obtained with the Skyrme energy-density functional, we propose an analytical 
expression for the description of the fusion probability, with which the measured evaporation 
residual cross sections can be reproduced acceptably well. Simultaneously, some special fusion 
reactions for synthesizing new elements 119 and 120 are studied. The predicted evaporation residual 
cross sections for ^"^Ti+^^^Bk are '--^ 10 — 150 fb at energies around the entrance-channel Coulomb 
barrier. For the fusion reactions synthesizing element 120 with projectiles ^^Cr and ^^Fe, the cross 
sections fall to a few femtobarns which seems beyond the limit of the available facilities. 
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Synthesis of super-heavy nuclei (SHN) through fusion reactions is a field of very intense 
studies in the recent decades [1-17]. Up to now, the superheavy elements Z = 107 ~ 118 
have already been synthesized [1-6] through "cold" fusion reactions that use lead or bismuth 
targets with appropriate projectiles following the emission of one or two neutrons from a 
"cold" compound system, or "hot" fusion reactions that use actinide targets from uranium 
to californium with beams of ^^Ca following the evaporation of 3 to 5 neutrons from a "hot" 
system. The reaction ^''Ti+^^^Cf for producing element 120 is currently being studied in GSI 
without a result up to now [18]. Besides the experiments, both the structure of SHN and the 
mechanism of fusion reactions are also being intensively investigated theoretically. On one 
hand, the precise calculation of the masses and the shell correction of SHN which plays a key 
role in determining the fission barrier and the center of the "island of stability" for SHN, is of 
great importance. Recently, Liu et al. [19] proposed an improved macroscopic-microscopic 
mass model (also called Weizsacker-Skyrme mass model) with an rms error of 336 keV with 
respect to the 2149 known masses and 248 keV to the measured a-decay energies of 46 
super-heavy nuclei, with which they found that the SHN with the largest shell corrections 
are located around Z = 116 ~ 120 and = 178 rather than at = 184. It is found that 
the difference of the calculated evaporation residual cross sections for reactions leading to 
element 120 reaches two orders of magnitude by adopting two different mass tables for SHN 
[15]. It is therefore interesting to investigate the production cross sections of SHN in fusion 
reactions by using the Weizsacker-Skyrme mass model. 

On the other hand, the fusion probability of heavy nuclei and super-heavy nuclei should be 
systematically investigated for testing the models and for giving reliable predictions of fusion 
reactions leading to new super-heavy elements. Theoretical support for these very time- 
consuming and extremely-expensive experiments is vital for choosing the optimum target- 
projectile-energy combinations and for the estimation of cross sections. In the practical 
calculation of the evaporation residue cross section, the reaction process leading to the 
synthesis of SHN can be divided into three steps. Firstly, the projectile is captured by the 
target and a dinuclear system is formed which then evolves into the compound nucleus, and 
finally, the compound nucleus loses its excitation energy mainly by the emission of particles 
and 7-rays and goes to its ground state. The simplified version of the evaporation residue 
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cross section is given by [14, 21] 

.m.) ''^cap 

Here, Ucap, -Pcn and Wsnr denote the capture cross section of the colhding nuclei overcoming 
the Coulomb barrier, the probability of the compound nucleus formation (i.e., the fusion 
probability) after the capture and the survival probability of the excited compound nucleus, 
respectively. The most unclear part may be the value of Pcn which is usually calculated 
by different models for the dynamics [10, 12, 13] based on the potential energy surface of 
the reaction system or by empirical formulas [12, 20]. In Refs. [21-24], the methods for the 
calculation of (Xcap and PVsur are well established in general, with the Skyrme energy- density 
functional and an empirical barrier distribution for describing the capture cross sections and 
with the HIVAP code [25-27] for describing the survival probability of the compound nuclei. 
Based on the previous work [21], the fusion probability in reactions leading to SHN will be 
further systematically investigated in this work by applying the Weizsacker-Skyrme mass 
model for describing the masses of SHN. 

We first investigate the relation between the fission barrier height, the quasi-fission barrier 
height and the evaporation residual cross sections for some fusion reactions leading to SHN. 
The fission barrier height B{ of a nucleus can be approximately estimated by the shell 
correction AE^^'^^^ and the macroscopic fission barrier Bf^'^'^ of the nucleus [21], i.e., 

Bi ^ 5f ^= + AE^^'^^ = + (Eexp - ^ld). (2) 

Here, Ecxp and E'ld denote the measured binding energy and the liquid drop energy of the 
nucleus (positive values), respectively. Fig. 1(a) shows the measured evaporation residual 
cross sections for some '^^Ca-induced reactions as a function of the charge number of the 
compound nucleus. The data are extracted from a figure in Oganessian's talk [28], and 
these data are also shown in Ref.[29]. At Z ^ 114—116 there exists a peak for the cross 
sections, which seems to be an evidence of the center of the "island of stability" for super- 
heavy nuclei being located at Z ^ 120. Fig. 1(b) shows the sum of the fission barrier 
height and the quasi-fission barrier height for some fusion reactions leading to the synthesis 
of SHN. Here, the quasi-fission barrier height B^f is defined as the depth of the capture 
pocket in the entrance-channel potential (see the sub- figure in Fig. 2), which is calculated 
by using the Skyrme energy- density functional [30] together with the extended Thomas- 
Fermi approximation including all terms up to the second order in the spatial derivatives 
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FIG. 1: (Color online) (a) Measured evaporation residual cross-sections of some ^'^Ca-induced 
reactions. The data are extracted from a figure in Oganessian's talk [28]. (b) Sum of fission 
barrier height and quasi-fission barrier height for some fusion reactions leading to the synthesis 
of super-heavy nuclei as a function of the charge number of SHN. The dot-dashed curves are to 
guide the eyes. The vertical dashed line shows the position oi Z = 118 and the squares on its right 
side denote the results for ^"^Ti-l-^^^Bk and ^'^Ti-l-^'^^Cf. (3) Quasi-fission barrier height for these 
reactions. The open and filled circles denote the results for "cold" and "hot" fusion reactions, 
respectively. 

(ETF2) as mentioned in [22]. For calculating the shell corrections in Fig. 1(b), the obtained 
binding energies and the corresponding liquid-drop energies -Eld of the SHN in Ref. [19] are 
adopted. One sees that the dependence of Bi + B^^ on the charge number is very similar 
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to the dependence of the evaporation residual cross sections on the charge number. The 
peak around Z ?^ 116 in Fig. 1(b) comes from the contributions of the shell correction and 
the quasi-fission barrier height. The quasi- fission barrier height decreases linearly with the 
charge number of the compound nucleus [see Fig. 1(c)], and the largest shell corrections 
are located around Z = 116 ~ 120 according to the Weizsacker-Skyrme mass model. The 
sum of the shell correction AE'^'^'^" and the quasi-fission barrier height results in the peak 
at Z 114 — 116 of the evaporation residual cross sections. Fig. 1 indicates that the quasi- 
fission barrier influences the formation of the super-heavy nuclei around the known "island of 
stability" in addition to the shell correction. For very asymmetric fusion reactions leading 
to intermediate and heavy nuclei rather than the SHN, the quasi-fission barrier height is 
relatively high. With the increase of the charge number of the projectile or target nuclei in 
the fusion reactions leading to the synthesis of SHN, the quasi-fission barrier height decreases 
gradually. 

It is known that the fission barrier height strongly influences the survival probability 
of super-heavy nuclei. In addition, one expects that the fusion probability increases with 
increasing of the quasi-fission barrier height and the incident energies in the reactions leading 
to the SHN [16]. To consider the contribution of the quasi-fission barrier, we propose an 
analytical formula for description of the fusion probability Pcn, 

Pcn{E*) = ^ exp (3fiqf + 0.3E*) . (3) 

Where, E* = Ec.m. + Q denotes the excitation energy of the compound nucleus (in MeV). 
In addition, we introduce a truncation for the value of Pcn? i-e. -Pcn should not be larger 
than 1. The parameters in Eq. (3) with C = exp(50|?7|) are determined by fitting the 
measured evaporation residual cross sections of ^^Ca-l-^'^^Cm and ^^Ca-|-^^^Cf. Here, t] = 
{Ai — A2) I (Ai -|- A'2) denotes the mass asymmetry of the reaction system. For ^*Ca-induced 
"hot" fusion reactions leading to SHN, \rj\ ^ 0.67. The calculation for the capture cross 
sections and the survival probability of compound nucleus are the same as those in Ref. [21] 
except that the masses of unknown nuclei are given by the Weizsacker-Skyrme mass model 
[19] rather than by the finite range droplet model (FRDM) [31]. Fig. 2 shows the fusion 
probability in reactions leading to SHN as a function of the quasi-fission barrier height 
Bq{. The solid squares denote the results for "hot" fusion reactions with an excitation 
energy of E* = 35 MeV. The open and filled circles denote the calculated results for "cold" 
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FIG. 2: (Color online) Fusion probability in reactions leading to super-heavy nuclei as a function of 
quasi- fission barrier height. The open circles and the filled circles denote the calculated results for 
"cold" fusion reactions with excitation energy E* = 15 MeV and the extracted fusion probability in 
Ref.[14] for the same reactions, respectively. The solid squares denote the results for "hot" fusion 
reactions with E* = 35 MeV. Sub-figure: Calculated entrance-channel potential for ^'^Ca-l-^'^'^Pb. 
-Bqf and Bq denote the depth of the pocket and the height of the Coulomb barrier, respectively. 

fusion reactions with E* = 15 MeV and the extracted fusion probability in Ref. [14] for 
the same reactions, respectively. Here the parameter C is slightly different for the "cold" 
fusion reactions which will be discussed later. The corresponding evaporation residual cross 
sections for these reactions will be shown in Fig. 3, Fig. 4 and Fig. 5. The shades show the 
estimated uncertainty of the calculated fusion probability with Eq.(3). With decreasing of 
the quasi-fission barrier height, the fusion probability decreases exponentially. The obtained 
fusion probabilities in this work and those in Ref. [14] are very close to each other for the 
"cold" fusion reactions. In addition, it seems that the obtained fusion probability for the 
"cold" fusion reaction and that for the "hot" fusion reaction are comparable when their 
quasi-fission barrier heights are close to each other. 
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FIG. 3: (Color online) Evaporation residual cross sections for some "hot" fusion reactions. The 
data are taken from Refs.[6, 32]. The dashed and solid lines denote the entrance-channel Coulomb 
barrier height Bq under the sudden approximation for the densities and the mean barrier height 
B^, respectively. The green and blue curves denote the corresponding calculation results for the 
3n and 4n channels with C = exp(50|r/|), respectively. The shades show the uncertainty of the 
calculated cjer for the 4n channel [except (h) in which the shades are for the 3n channel]. 

Fig. 3 shows the evaporation residual cross sections for some "hot" fusion reactions. The 
open triangles and solid circles denote the experimental data for the 3n and 4n channels, 
respectively. The green and blue curves denote the corresponding calculated results. The 
dashed and solid lines denote the entrance-channel Coulomb barrier height Bq from the 
Skyrme energy-density functional and the mean barrier height ^ 0.956-Bo according 
to the empirical barrier distribution [22], respectively. The shades show the corresponding 
uncertainty of this approach due to the systematic errors in the calculation of Ccap, W^sur 
and PcN, which is about a factor of 4.4 for o"er at energies above the mean barrier height 
i?ni- The estimated uncertainties [21] for (Xcap, W^sur and Pcn are about factors of 1.18, 1.85 
and 2.0 at -Ec.m. > B„„ respectively. We would like to state that the parameters in Eq.(3) 
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FIG. 4: (Color online) The same as Fig. 3, but for some "cold" fusion reactions. Here, we adopt 
C = exp(47|7y|) in the calculations. The data are taken from Ref.[32]. The shades show the 
uncertainty of the calculated cjer for the In channel. 

are determined just by reactions ^®Ca+^^^Cm and ^'^Ca+^^^Cf. With the same parameters 
we find that the experimental data of other "hot" fusion reactions can also be reproduced 
reasonably well. 

To further test the model, the "cold" fusion reactions leading to SHN are also investigated 
systematically. In Fig. 4, we show the evaporation residual cross sections for some "cold" 
fusion reactions. The solid squares and open circles denote the experimental data for the 
In and 2n channels, respectively. The blue and green solid curves denote the corresponding 
calculated results. Here, we slightly change the parameter C = exp(47|?7|) considering the 
differences between the "cold " fusion reactions and the "hot" fusion reactions. Carrying 
out the same calculations as for the "hot" fusion reactions, we can roughly reproduce the 
measured evaporation residual cross sections for the "cold" fusion reactions. The larger 
uncertainty for the "cold" fusion reactions at sub-barrier energies is due to the factor g in 
the barrier distribution function for calculating the capture cross sections [21]. 
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FIG. 5: (Color online) The same as Fig. 3 but for fusion reactions leading to the synthesis of 
elements 119 and 120 (in femtobarn). The green solid and blue dot-dashed curves denote the 
predicted results for the 3n and 4n channels with C = exp(50|7/|), respectively. 

TABLE I: The predicted optimal evaporation residual cross section (in femtobarn) for four reac- 
tions leading to elements 119 and 120. 



Reference 


50xi+249Bk 


50Ti+249Cf 


54Cr+248Cm 


58pg^244p^j 


Liu [10] 


600 


~ 100 






Nan Wang [34] 


~ 1000 


~ 200 


~ 40 


~ 30 


Zagrebaev [12] 


~ 50 


~ 40 


~ 20 


~ 5 


This work 


- 35 


~ 20 


~ 5 


~ 3 



The above calculations give us great confidence for investigating the evaporation residual 
cross sections of fusion reactions leading to the synthesis of new elements. Fig. 5 shows 
the predicted evaporation residual cross sections for the reactions a^Ti+^gyBk, g^Ti+^glCf, 
24Cr-|-^9gCm and 26Fe+^g|Pu leading to the synthesis of the elements 119 and 120. The 
solid and dashed curves denote the results for the 3n and 4n channels, respectively. For the 
reaction ^°Ti -|-^^^ Bk — 119-1- 4n, the obtained evaporation residual cross sections o"er 



9 



are about 10 ~ 150 femtobarn at incident energies around the entrance-channel Coulomb 
barrier Bq. For the reactions synthesizing element 120, the calculated cter are smaller than 
those for element 119, and the o"er falls to a few femtobarns for the latter two reactions 
^^Cr+^^^Cm and ^^Fe+^'^^Pu. Here, we would like to state that the parameter C in Eq.(3) 
for the description of the fusion probability depends on the parameter sets for calculating 
the survival probability of super-heavy nuclei. If one takes a different liquid-drop formula 
for calculating the fission barrier, the parameter C must be re-adjusted. Adopting different 
formulas for describing the liquid-drop energies [19, 33] and the binding energies of yet 
unmeasured nuclei [19, 31] (the parameter C is re-adjusted to fit the measured o"er of "hot" 
fusion reactions), we find that the calculated evaporation residual cross sections for the four 
reactions leading to the elements 119 and 120 are close to the corresponding results in Fig. 
5, with deviations which are smaller than the uncertainty of this approach. To compare 
with the results from other models, we list in Table I the calculated optimal evaporation 
residual cross sections from four different models for the four reactions mentioned in Fig. 
5. One sees that for ^°Ti-|-^^^Bk and ^°Ti-|-^^^Cf the results from Ref.[10] with the fusion- 
by-diffusion model and Ref.[34] with the dinuclear system model are in the same order of 
magnitude (hundreds of fb), while the results from Ref.[12] and those in this work are in 
the same order of magnitude (tens of fb). All models in Table I predict that the optimal 
evaporation residual cross sections for ^°Ti -|-^^^ Cf — > 299-a;n ]^20 are smaller than those for 

50Ti+249 299-xn ^^g^ 

In summary, the entrance-channel Coulomb barrier determines the capture cross sections, 
the quasi-fission barrier due to the pocket in the entrance-channel potential influences the 
fusion probability and the fission barrier is related to the survival probability, all play an 
important role for a reliable calculation of the evaporation residual cross sections in fusion 
reactions leading to the synthesis of SHN. Of course the dynamics of the fusion process and 
the zero-point motion for fission process [35, 36] also influence the cross sections significantly, 
but it is beyond the scope of this work. Based on the previously established methods for 
the calculation of capture cross sections and the survival probability, we further proposed 
an analytical formula with only three fixed parameters for a systematic description of the 
fusion probability in reactions leading to SHN. With this formula the measured evaporation 
residual cross sections for the "hot" fusion reactions can be reproduced reasonably well. 
The predicted evaporation residual cross sections cter for the reaction ^°Ti-|-^^^Bk are about 
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10 ~ 150 femtobarn at incident energies around the entrance-channel Coulomb barrier, and 
the predicted cter for the reactions ^^Cr+^^®Cm and ^'^Fe+^^^Pu falls to a few femtobarns 
which seems beyond the limit of the available facilities. 
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